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Abstract 

We study the Boltzmann-Langevin equation which describes the dynamics of hot Yang- 
Mills fields with typical momenta of order of the magnetic screening scale g 2 T. It 
is transformed into a path integral and Feynman rules are obtained. We find that 
the leading log Langevin equation can be systematically improved in a well behaved 
expansion in log(l/g) -1 . The result by Arnold and Yaffe that the leading log Langevin 
equation is still valid at next-to-leading-log order is confirmed. We also confirm their 
result for the next-to- leading-log damping coefficient, or color conductivity, which is 
shown to be gauge fixing independent for a certain class of gauges. The frequency scale 
g 2 T does not contribute to this result, but it does contribute, by power counting, to the 
transverse gauge field propagator. Going beyond a perturbative expansion we find 1- 
loop ultraviolet divergences which cannot be removed by renormalizing the parameters 
in the Boltzmann-Langevin equation. 
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1 Introduction 



We consider hot non-abelian gauge fields at a temperature T sufficiently large so that 
the the running coupling g = g(T) is small. Even then gauge fields with typical mo- 
menta of order g 2 T, the so called magnetic screening scale, are strongly coupled fflj. 
Despite this, physical quantities have an expansion in powers of g times possible loga- 
rithms of g. Only the coefficients in this expansion cannot be calculated perturbatively 
when they receive contributions from magnetic scale momenta. For example, the or- 
der g 6 contribution to the free energy comes with a non-perturbative coefficient. This 
coefficient can be calculated by a lattice simulation of 3-dimensional pure Yang-Mills 
theory, which was shown by using dimensional reduction [Q]. This technique can only 
be applied to static quantities, like the free energy or equal time correlation functions. 

In this paper we consider dynamical quantities, which are determined by unequal 
(real) time correlation functions. While the free energy is sensitive to the magnetic 
scale only at relatively high order in the weak coupling expansion, there are dynamical 
quantities which are determined by this scale at leading order. An important example is 
the Chern-Simons diffusion rate in unbroken non-abelian gauge theories, often referred 
to as the hot sphaleron rate. In the standard electroweak theory it determines the 
rate for anomalous baryon number violation above the electroweak phase transition or 
cross-over temperature T c ~ 100 GeV It is a crucial ingredient in scenarios which try 
to explain the baryon asymmetry of the universe, e.g., through leptogenesis || or elec- 
troweak baryogenesis [ffl. In QCD Chern-Simons diffusion leads to non-conservation 
of quark helicity even in the chiral limit |J. Another physical quantity which is sensi- 
tive to the magnetic screening scale at leading order is the friction experienced by an 
electroweak phase transition bubble wall ||. 

Euclidean lattice simulations are of limited use for calculating real time observables (], 
and real time simulations of quantum field theory are not possible. However, even for 
dynamical quantities one can use perturbation theory to integrate out short distance 
modes to obtain an effective theory for the non-perturbative long-distance dynamics. 
The first step is to integrate out "hard" physics associated with virtual momenta of 
order T. At leading order in the gradient expansion one obtains the so called hard 

2 Close to T c ~ 100 one must take the Higgs field into account unless the transition is strongly 
first order. 

3 Recently methods have been developed for obtaining real time correlation functions from Eu- 
clidean lattice simulations using the so called maximal entropy method However, this cannot be 
applied to Chern-Simons diffusion which is exponentially suppressed in Euclidean time. 
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thermal loop effective theory || . The remaining degrees of freedom are gauge fields with 
characteristic wave vectors k< gT. Their occupation number n(|k|) = l/(e' k '/ T — 1) ~ 
T/|k| is large which means that they can be treated as classical fields. This opens 
the possibility for a non-perturbative treatment because classical field theories can be 
treated on the lattice in real time fllQfl . The hard thermal loop effective theory was 



used for non-perturbative lattice simulations (TT|-[l3j which, however, do not have a 
continuum limit fl4}| . 

One can, however, go further using perturbation theory. Integrating out all frequency 
scales and spatial momenta f] larger than g 2 T one finds that at leading order in g the 



dynamics of gauge fields is described by the Boltzmann-Langevin equation [15 



(C + v D)W + -4-(ci + v • D)v DxB = vE + «£, (1.1) 



which is a classical field equation for A and W. D^ b = <5 a fc<9 M — gf a b c ^-c * s the covariant 
derivative in the adjoint representation, E l = F l ° and B % = —±e l i k Fi k are the non- 
abelian electric and magnetic fields, mo is the Debye or electric screening mass which 
is of order gT, and d denotes the number of spatial dimensions. We choose d = 
3 — e to regularize ultraviolet and infrared divergences. W(x, v) transforms under the 
adjoint representation and represents color charge fluctuations of plasma particles with 
momenta of order T and velocity v. We neglect fermion masses. Then these particles 
move at the speed of light so that v 2 = 1. The 4-vector v is defined as = (l,v). 
C is a linear collision operator of order g 2 T and depends logarithmically on the cutoff 
separating the scales gT and g 2 T. It is local in space and acts only on the velocity 
variable, 

C/(v)= / C(v,v')/(v'). (1.2) 

Jv' 

Here J* (• • •) = Q^ 1 J d d ~ 1 v(- ■ •) is the normalized integral over the d — 1 dimensional 
unit sphere, Q = An for d — 3. C commutes with rotations of v. Therefore it can be 
diagonalized by expanding W(x, v) in spherical harmonics f] YJ m (v), and the eigenvalues 
depend only on I. The I = eigenvalue vanishes, and the coefficient c\ in Eq. ( |1 . 1| ) is 
the I = 1 eigenvalue. W contains only spherical harmonics with respect v with / > 2, 



4 Since Lorentz invariance is lost at finite temperature, frequencies and spatial momenta have to be 
treated separately. 

5 Or their generalization in d ^ 3 dimensions. To keep the terminology simple, we refer to them as 
spherical harmonics. 



2 



i.e., 



v"W(x,v) = 0. (1.3) 



£(x, v) is a Gaussian white noise with the 2-point function 

9T 

(Ux,v)Ux',v')) = — C (v,V)5 ab 5 d+1 (x-x'). (1.4) 
m D 

A non-trivial check of Eq. ( |1 . 1| ) is to see whether it reproduces the known thermo- 
dynamics of the magnetic scale gauge fields. In Ref. [TJ|] it was indeed shown that 
Eq. ( Jl.l| ) generates a thermal ensemble with a Boltzmann weight e~ H ^ T where 

H = ~ J d 3 xi^B 2 + ml y^ 2 }- ( L5 ) 

Thus for equal time correlation functions W acts like a free field and does not affect 
the gauge fields. The term |B 2 is just the Lagrangian of the dimensionally reduced 
theory for the scale g 2 T [|J. 

In the leading log approximation, neglecting terms suppressed by powers of [log(l/g)] _1 , 
Eq. ( [LTD can ^ e replaced by a much simpler equation of motion which contains only 
the gauge fields fll6| , 



D x B = 7E+ C (1.6) 

Here £ is again a Gaussian white noise which is proportional to the I = 1 projection 
of £. The damping coefficient 7, the color conductivity, is proportional to T / log(l/g). 
The only spatial momentum scale in Eq. ( |1.6|) is k ~ g 2 T. Estimating the LHS as k 2 A 
and the first term on the RHS as jk°A one can see that the only frequency scale in 
Eq. (|1.6j ) is k ~ g 4 T (modulo logarithms of the coupling). The Boltzmann-Langevin 
equation (|1 . 1|) , on the other hand, contains two frequency scales, g 4 T and g 2 T, again 
modulo logarithms [15 j. The question whether the frequency scale g 2 T affects the 
non-perturbative gauge field dynamics has not been addressed so far. 

Eq. ( |l.b1 ) solves the notorious ultraviolet problems mentioned above. It is ultravi- 



olet finite and has been used to calculate the hot sphaleron rate by Moore ||17|| . Its 
obvious disadvantage is that it's only valid in the logarithmic approximation. The 
Boltzmann equation which was obtained by integrating out the momentum scale gT 



1(| f\ and which served as a starting point for obtaining Eqs. (|1.6|) and (|1 . 1|) is valid to 



For a Feynman-diagram interpretation, see [18l [L9| 
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leading order in g and all orders in [log(l/g)] 1 , but it has the same 1-loop ultraviolet 



divergences as classical Yang-Mills theory. One can hope that Eq. (|1 . 1|) has a better 
ultraviolet behavior because it does not contain propagating gauge field modes which 
are responsible for these divergences. 

An equation similar to ( |1 . 1| ) has been obtained by Arnold P0 |. The difference be- 
tween Arnold's Langevin equation and Eq. ( |1 . 1| ) is that the former does not contain 
a time derivative of W. Thus in |20| W is not dynamical, but it is fixed by A and 
£ at the same instant of time. The only dynamical degrees of freedom are the gauge 
fields. On first sight this difference appears to be minor. However, if one discretizes 
time for solving these equations on the lattice with temporal lattice spacing Ax° the 
two equations require rather different treatments. For Eq. ( |1.1| ) A and W at time 
x° + Aa; and spatial coordinate x are determined by A, W, a finite number of their 
derivatives and £, all taken at x M = (x°,x). This procedure is exactly the same as the 
one one would use for solving Eq. ( |1.6| ) ]T7|] . For Arnold's Langevin equation, on the 
other hand, one first has to solve for W at time a; , so that W(x°) becomes a spatially 
non-local functional of A and £ at time x°. Then one can determine A(x° + Ax°) from 
the so obtained W(x°) and from A(x°) and £(x ). Since W is linear in £, Arnold's 
Langevin equation can be viewed as a spatially non-local generalization of Eq. ( |1.6| ) in 
which the color conductivity 7 and the noise £ are gauge field dependent. The time 



discretization of such a Langevin equation is ambiguous, but it was argued in Ref. [2C 
that the ambiguity can be fixed by demanding that the equation yields the correct 
thermodynamics of magnetic scale gauge fields. In the path integral formulation of 
this equation (cf. Sec. 0) it is necessary to introduce additional terms in the action []. 
No such complications arise for Eq. ( |1 . 1| ) . 

Arnold and Yaffe used the Langevin equation of [20| to show that Eq. ( |1.6| ) is still 
valid at next-to- leading logarithmic order (NLLO), if one uses the next-to-leading log- 
arithmic order (NLLO) color conductivity 7, which was calculated in Ref. [21|. It was 
argued that the additional terms in the path integral action mentioned above do not 
contribute at next-to-leading logarithmic order. 

The purpose of this paper is three-fold. First, we want to study the large momentum 
behavior of Eq. ( |1 . 1| ) to see whether it is renormalizable. This is necessary if one wants 
to use Eq. (LI) for non-perturbative lattice simulations and take the continuum limit. 
Secondly, we would like to understand whether Eq. ( |1.1| ) gives the same gauge field 
dynamics as the Langevin equation obtained by Arnold ||20|| . Finally we calculate the 

7 These terms are rather complicated. To quote the author of Ref. M, 'they are ugly as sin'. 
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color conductivity at next-to-leading logarithmic order. We confirm the result obtained 
by Arnold and Yaffe [21], and we investigate its gauge fixing dependence. 

To address these issues it is convenient to write the equation of motion ( |1.1| ) as a 
path integral. This method is briefly reviewed in Sec. |2] for a non-gauge theory. To 
apply it to Eq. ( |1.1| ) one has to fix the gauge which is discussed in Sec. [| The path 
integral representation and the corresponding Feynman rules for Eq. ( |1.1| ) are obtained 
in Sec. In Sec. |5] we study the large momentum behavior of the various propagators. 
We shall find 1-loop ultraviolet divergences which cannot be removed by renormalizing 
the parameters of the action (Sec. |5.1| ). In Sec. ^| we discuss the logarithmic approxi- 
mation beyond leading order. Following Arnold and Yaffe we consider Wilson loops in 
Sec. [7| at next-to-leading logarithmic order (NLLO) to compute the NLLO color con- 
ductivity. In Sec. H we study the effect of the modes with k° ~ k ~ g 2 T on the spatial 
gauge field propagator. Sec. [| contains a summary and a discussion. In Appendix |A] we 
list some properties of the operator G defined in Eq. ( 4.15 ). Appendix [B| describes the 
power counting for propagators and vertices used in Sects. [7| and ||. It also contains the 
approximated expressions for the propagators which are used in Sec. 0. The frequency 
integrals needed for the 1-loop selfenergies are listed in Appendix |C|. Appendix |B] con- 
tains the results for the individual 1-loop diagrams contributing to the A propagator 
which gives the dominant contribution to the Wilson loops discussed in Sec. [7]. 
Notation 4-vectors and 3-vectors are denoted by italics and boldface, respectively. 
The signature of the metric is ( H ). 



2 Path integral representation of Langevin equa- 
tions 

To investigate the properties of a Langevin equation it is convenient to re-write it as 
a path integral [p^l , which allows the use of the standard methods of quantum field 
theory. In this section we illustrate this method for a non-gauge theory. The fields are 
denoted by (p l and i represents all indices including space and time variables. Similarly, 
we collectively label the components of the equation of motion with the index a, 

S a [ip]=C- (2-1) 
We write the 2-point function of the noise as 

(£ a £P) = 2r5 a/3 . (2.2) 
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Physical observables are obtained by averaging over the noise with a Gaussian weight 
(O) = f [d£] exp j-^CC } 0(<p sMtion ), (2.3) 



where (/'solution is the solution to the equation of motion (|2.1|). Now one introduces an 
integral over if via 



1 = J [ dv ]det(j^j5(£-d. (2.4) 

The 5-function in ( |2.4| ) is represented by an integral over a Lagrange multiplier field A, 

(O) = J [cW£dA]det (^pj exp {-^FT + *A Q (£ a ~ £ Q )} 0[<p). (2.5) 

Finally one integrates over the noise and introduces ghost and anti-ghost fields rj and 
fj to represent the determinant which gives 



(0)=l[*pd\dndfj]e->0[<p]. (2.6) 

The action 

S = r\ a \ a - i\ a £ a - ^—rf (2.7) 
is invariant under the BRST-transformation 

V = erf, 5f] a = — isX 01 , 5X a = 8rf = (2.8) 

with the ant i- commuting parameter e. This procedure can be applied to gauge theories, 
but first one has to fix the gauge. 

3 Gauge fixing 

Gauge covariant Langevin equations determine the gauge fields only up to a gauge 
transformation. Thus in order to use Eq. (|2.4| ) to introduce a path integral one first 



has to fix a gauge. Zinn- Justin and Zwanziger |23j showed that the gauge condition 

A°(t) = v[A(t)} (3.1) 
does not affect expectation values of gauge invariant operators 0[A]. A small variation 

v[A] -> v[A] + e[A] (3.2) 
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can be compensated by a gauge transformation of A. It is easy to see that this result 
also holds for gauge invariant operators which depend both on the spatial gauge fields 
A and on A . 

A special case of Q3.1| ) is 

kA° + V ■ A = 0, (3.3) 

which smoothly interpolates between Coulomb (« = 0) and temporal gauge (k — > oo). 
The free fields in this gauge are related to A in Coulomb gauge by 

k 2 i 

A ° {k) = TFTikV^^^' 

t-A(fc) = fc0 J^2 //t A °Coui(fe) (free fields). (3.4) 
With this relation we can obtain the propagators for the gauge (|3.3|) from the Coulomb 



gauge propagators listed in Sec. 4.1 



To obtain a path integral representation one can simply use the gauge condition 
( |3.3|) to eliminate A from the equations of motion, Eq. ( |1 . 1| ) or Eq. (|1.6|) , and then 
apply the procedure of Sec. ||. Closed ghost loops vanish in this case, so they do not 
contribute to amplitudes of the fields (p and A. 

A slightly different method was used by Arnold and Yaffe [21] who used temporal 
A = gauge to obtain the path integral representation. For the leading log Langevin 
equation ( |1.6|) this leads to a supersymmetric action in the path integral. It is not 
clear, however, whether the analogous action for Eq. (|1.1|) would be supersymmetric. 
Then the Fadeev-Popov procedure was used in |H] to switch to Coulomb gauge. 



4 Path integral representation of the Boltzmann- 
Langevin equation 

We can now write the path integral representation of Eq. ( |1 . 1| ) as 

J[dAdWd\dr]df]}e- s (4.1) 



with the action 

' 1 ' (C + vD)W 



S = I d*X I i -ZrXCX - iX 



'v I ml 



+ -4-(ci + v ■ D)v D x B - v E 



ml 



+ Sghost- (4.2) 
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As was mentioned above, closed ghost loops vanish in this theory. In this paper we are 
only interested in amplitudes of the physical fields A and W and of A. Therefore we 
can ignore the ghost part of the action Sghost m t ne following. 
The vertices can be read off from the interaction Lagrangian 



Oil 



-igf abc I (2n) d+1 6(k 1 + k 2 + k 3 ) 

' fcl,fe2,fc3 



X a (kt, v) [A° b {k 2 ) - v • A b (k 2 )]W c (k 3 , v) 



1 d 

+ 



[-ic x \\{k x )-W x Xl{k x )\ 
x \{k 2 - k x ) n 5 u + (k 3 - k 2 y5 nl + (k x - h) l 5 in ] 



+ Af {h) [6 nm P: l (k 2 )k 2 2 - 6 lm Pr(k 3 )kl] 



A a (fc 1 )-A fe (A; 2 )A0(A;3) + O(AA 3 ) 



A l b (k 2 )A n c (k 3 



(4.3) 



with 



d d+l k 



Furthermore, we have introduced the 3— vector A through 



A^ = / u*A(v) 



and the traceless tensor 



X ij = 



v l v 3 - -5 lJ ) A(v), 



(4.4) 



(4.5) 



(4.6) 



which represent the / = 1 and / = 2 projections of A(v), respectively. Finally, the 
transverse projector is 

Pl\k) = 5 ij - titi . (4.7) 

In Eq. ( [4.3[ ) we have not displayed the AA 3 and AA 4 vertices which do not play a role 
in the following. The only interaction in Eq. ( f4.3| ) which contains derivatives is the 
AA 2 vertex. For power counting we can estimate this vertex as 

cJkl 



(AAA) ~ g- 



2 ' 
D 



(4.8) 
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where k is a typical spatial momentum carried by one of the fields. The remaining 
non-derivative interaction vertices are simply of order g. 



4.1 Propagators 



Obtaining the propagators for Eq. ( |4.2| ) is tedious but straightforward. To simplify the 
notation it is useful to think of C and v as operators acting on the space of functions 
on the 2-sphere. These functions are represented by bras and kets. We define |v) as 
the eigenvector of the operators v l (i — 1, . . . , d) with eigenvalues v\ 



and with the normalization 



v l \v) = v % \v) 



(vIv'HO^V-v') 



(4.9) 



(4.10) 



Here S^ s ^ is the delta-function on the d— 1 dimensional unit sphere, the surface area 
of which is denoted by Q. Then the unit operator is 



It is also convenient to introduce 



These are normalized such that 



|v)(v| 



) = / ^1 V >- 



( u | u ) = l, = S 13 . 



d 



(4.n; 



(4.12) 



(4.13) 



The action (|4.2 ) can then be written as 



S 



d A xl -^-(A|(7|A) 



m 



D 



(C + v-D)\W) 



m 



D 



(4.14) 



To determine the momentum space propagators one has to invert C — iv-k. We denote 
its inverse by 

-l 



G(k) = 



C — iv ■ k 



(4.15) 
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For writing propagators involving gauge fields it is convenient to introduce 



which is symmetric under exchange of /i and v. Its spatial components are decomposed 
into transverse and longitudinal parts, 

G ij = Pi j G t + k r k j G e . (4.17) 

Some properties of G and related functions are listed in Appendix A. 

We give the expressions for propagators in Coulomb gauge. The propagators for 
general flow gauge can then be obtained by using Eq. Q3.4Q . We denote spatial gauge 
fields by wavy, A by dashed, W by full, and the Lagrange multiplier fields A by dotted 
lines. 

The AA propagator vanishes, 

A. • A = 0. (4.18) 

This is an immediate consequence of the invariance under the BRST transformation 
( |2.8D . Due to ghost-number conservation we have (At/) = 0. Performing a BRST 
transformation in the path integral for (At/) = then gives (AA) = 0. 

The off-diagonal propagators mixing the physical fields with A are not invariant 
under changing sign of the momentum. Therefore we explicitly display the momentum 
variables, 

A%k)m • \(-k,v) = i(S*{k)\V), (4.19) 



A°{k)+ .A(-ife,v) = z<5°(A;)|v). (4.20) 

The expressions 

(S\k)\=mlA(k)Pi\k)(>\G(k), (4.21) 

with 

A(Jfe) = [1 + idkoGtik)}" 1 A(Jfe), (4.22) 



A(*) = 



k 2 - tmlk G t (k)[l + idk G t {k)]' 



(4.23) 
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and 



<S°(*)| 



k*C\G(k)(l-P ) 



(4.24) 



will also appear in the A^ - — W propagators below. The W — A propagator can be 
written as 

W(k,v)m .A(-£;,v')) = i(v\(l-P -Pi)s(k)(l-P^\V). (4.25) 

P and Pi are the projectors onto the / = and I = 1 sectors, 

p = rx°i, A = ^r)(i 

and 



(4.26) 



S(k) = G(k)- 
1 



-G^l^tftf ( j \G(k) 



G t (k) 



G e (k) 

l-k 2 A(k) G(k)\*)P?(k){>\G(k). 



(4.27) 



Note that all propagators which mix A with the " matter" fields and W have poles 
only in the lower half of the complex k° plane. As a consequence, loops which contain 
only such propagators vanish. 

To obtain the propagators for the physical fields, one uses 



G(k)CG(—k) = - G(k) + G(-k)\, 



and 



[G t (k)+G t (-k)]A(k)A(-k) 



mfy k 



[A(fc) - A(-fc)] 



After some algebra one finds that the gauge field propagators take the form 



A 1 



Ai = TP: j -A(k) + (k^-k), 



A°+- — .A° = 



T 1 

<k 2 G,(fc) 



(4.28) 
(4.29) 

(4.30) 
(4.31) 



Note that 



A*(fc) = A(-fc). 



(4.32) 



11 



Figure 1: Examples for 1-loop diagrams which are ultraviolet divergent, and for which 
no counterterms are available in the action (|4.2j ). 



The off-diagonal propagators are 

W • wA* = -^r(v|(l- P 1 )\S\k)) + (k ~^ -k), (4.33) 

mi 

W. • A = _L( v |(i_A)|S°(fc)> + (*;-►-*;), (4.34) 

where \S^) has been defined in Eqs. ( |4.21| ), ( |4.24| ) above. Finally, using again Eqs. ( |4.28| ), 
( [4.29| ) one finds the W propagator as 

T 



W m • W 



ml 



(v| (l - Po - A) S(k) (l-Po- A) |V) + (k -k). (4.35) 



5 Large momentum behavior of the propagators, 
1— loop ultraviolet divergences, and renormali- 
zability 

In this section we first study the large momentum behavior of the propagators to see 
what kind of ultraviolet divergences can appear. The large momentum behavior of the 
A propagator, Eq. ( f4.31p , can be easily estimated. We have Gg ~ 1/k, so that 



(A°A°) ~ i. (5.1) 



The W propagator (|4.35| ) shows the same behavior, 



(WW) ~ -. (5.2) 
k 



Now consider the transverse gauge field propagator (|4. 30|) 



(AA) ~ 1 [A(*) - A(-fc)] . (5.3) 
12 



One can write the square bracket as 

A(k)-A(-k) = 2zm|fe°Regt fg 4) 

[k 2 + /^(m 2 , - dk 2 )lmG t ] 2 + fc 2 (m 2 - dk 2 )ReG t 2 ' 

For k 2 ^> one can simplify the denominator, 

Am A/ 1 2imlk°ReG t 

A(k) - A(-k) ~ —r . (5.5) 

k 4 [1 - ^ ImG t ] 2 + c^ReG 2 

For large momenta the products k°KeGt and /c°ImG t are at most of order unity. Thus 
for k 2 ^> m 2 -, the transverse gauge field propagator (|4.30| ) scales like 

(AA) ~ ^4- (5.6) 

It is a little puzzling that this asymptotic behavior sets in only when k 2 S> mf,, which 
is well outside the range of validity of Eq. ( |1.1|) , k 2 <C m^. 

A consequence of Eq. ( |5.6|) is that the diagrams (p.2| ), (|D.3|) , ( p.5|) in Appendix 0, 



which contain transverse gauge field propagators, would be ultraviolet finite when 
computed with the propagator ( |4.30|) . The results in Appendix [D] are obtained by 
using the approximated propagators of Appendix || giving ultraviolet divergent results 
(see the discussion at the end of Sec. ^j. 

5.1 Renormalizability 

Now we address the question whether Eq. ( |i.lp is a renormalizable effective theory. 
As we have seen above, certain propagators, like (WW) and (A A ), fall off only like 
1/k for k —>■ oo. Therefore diagrams with 3 and 4 external A lines, which contain only 
these propagators (see Fig. |5|), are linearly and logarithmically divergent. The terms 
in the action ( [4.2[ ), on the other hand, are at most quadratic in A. Thus in order to 
renormalize these divergences one would have to introduce new terms in the action. 
But then the path integral ( [4.1| ) would no longer be equivalent to the Langevin equation 
Thus it would not be possible to use the theory ( |1.1| ) for lattice simulations if 
one wants to take the continuum limit If one is interested in calculating corrections 
to the leading log Langevin equation in an expansion in [log(l/g)] -1 this does not pose 
a problem since the divergent terms would contribute at higher orders in g. 

8 It is of course possible that these divergent contributions cancel among different diagrams. How- 
ever, I do not see a reason why such a cancellation should occur. 



13 



6 The logarithmic approximation beyond leading 
order 



Arnold and Yaffe have argued that the leading log Langevin equation (|1.6j ) is also valid 
at next-to- leading logarithmic order (NLLO). We will now present a simple argument 
why this should indeed be the case. 

Let us first recall the argument which takes us from Eq. (|1.1|) to Eq. ( P..6Q in the 
leading log approximation. 

For the non-perturbative fields with spatial momenta of order g 2 T the collision term, 
which is of order g 2 log(l/g)T, is large compared to the term containing v • D. If one 
neglects the latter, the I = 1 projection of Eq. ( |1 . 1| ) does not contain W and is thus 
a closed equation for the gauge fields alone 0, which is precisely Eq. ( |1.6|) . Similarly, 
the equations for different I and m components of W decouple in the leading log 
approximation. 

Now imagine continuing this approximation scheme by expanding in powers of L _1 = 
log(l/g) _1 . For the moment we ignore the Fourier components with k ~ Lg 2 T. At 
NLL order the / = 1 projection of Eq. ( |1.1[ ) contains a term which schematically looks 
like v • ~DWi = 2- At this order we can use the leading log approximation for W\ = i to 
write 

v ■ BW l=2 ~ v • DCr 1 ^ ~ logil/g)' 1 ^. (6.1) 

Thus it would appear that the correction term is suppressed by only a single power of 
L^ 1 . However, since the correlation functions of £ are diagonal in the (l,m) basis (see 
Eq. ( |1.4| )), one needs at least two factors of £;=2 to obtain a non- vanishing contribution. 
But then one also gets two factors of L^ 1 , which shows that the correction term in 
Eq. (|6.1|) would become relevant only at next-to-next-to-leading log order. Thus even 



at NLLO we can ignore the the I > 2 projections of Eq. (|4.2|) , and we can drop W field 
from fl4.2|) . Then we obtain the action 

S h = J d 4 x j-^ A 2 - i X • x B - | (6.2) 

with 

7b = (6.3) 
dc\ 



3 For a detailed discussion see fE4l 
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The modes with spatial momenta of order of the collision term, k ~ g 2 LT, Q enter 
the soft non-perturbative dynamics through loops. They can be treated perturbatively, 
where the loop expansion parameter is They change the coefficients in the action 
which becomes 

S = J d 4 x^Z x A 2 - % X (Jj-Z B D x B - Z E E^j J . (6.4) 



The k ~ C modes will also introduce higher dimensional operators in Eq. QS.4p . These, 
however, are suppressed by powers of g 2 T/C ~ L" 1 . 

The Lagrange multiplier field A can be rescaled to bring the action back into the 
form O, 

A -> Z E X A, (6.5) 

so that Eq. (|6.4| ) becomes 

S = [ d 4 x {-^ X 2 - t X - ( -^B xB - e] \ . (6.6) 

Now we can read off the color conductivity 7 which takes into account of the k ~ g 2 LT 
modes from the coefficient of the second term in Eq. ( |6.4|) , 



7 = f^7b- (6.7) 



Alternatively one could read off 7 from the first term in Eq. (|6.4|) because we already 
know that both Eq. fl6.4Q and Eq. (|1 . 1| ) generate a thermal ensemble of magnetic scale 
gauge fields at temperature T [^2[ ^] (see also Eq. (|1.5|) ). In other words, the 
temperature T in Eq. ( |6.6| ) is the same temperature appearing in Eq. ( |i.lp . Therefore 
we must have 

7 = ^7b- (6.8) 

Thus the Z-factors must be related by 

Z x = Z E Z B , (6.9) 

and the action contains only one unknown physical parameter 7. Thus in order to 
determine 7 rather than to calculate the Z-factors it is sufficient to calculate one single 
quantity in the theories (|1.1| ) and ( |1.6| ) and require that the results match. 
10 For simplicity we will oftentimes write k ~ C instead of k ~ g 2 LT. 
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We have seen in Sec. ^| that loops containing gauge field propagators in general 
receive non-trivial contributions from momenta of order m-Q ~ gT. These will not be 
included in the calculation of 7 in Sec. 0, i.e., the calculation in Sec. is performed with 
an ultraviolet cutoff separating the scales gT and g 2 T. Technically this is achieved in 
dimensional regularization by not using the propagators of Sec. |]but the approximated 
expressions for k 2 <C m^. These are listed in Appendix [FJ. 



7 Determining 7 by calculating Wilson loops 

To determine the next-to-leading log color conductivity 7 we follow Arnold and Yaffe 
[ PH and match the results for a rectangular Wilson loop W(t, R), with spatial extent R 
and temporal extent t, in the theory ( |6.4| ) and in the underlying "microscopic" theory, 
which in our case is Eq. ( |1 . 1[ ) . The Wilson loop is defined as 

W(t,R) = djHr^Pexp ji^ J dx ■ AjJ , (7.1) 

where P denotes path ordering, and the integration is along a closed contour starting 
at x M = 0, then going a distance R in the z- direction, then a distance t in the time 
direction, back in the ^-direction to z = and finally back to x M = along the time 
direction. The normalization factor d^ 1 , where g?r is the dimension of the representation 
R associated with the Wilson loop, is included so that \W\ < 1. 

One has to choose R such that R^ 1 3> g 2 T, so that one can obtain perturbative 
contributions to W (see Sec. |7.1|). Furthermore, one needs R^ 1 g 2 log(l/g)T, where 



both theories ( p..l|) and ( |1.6|) are valid. Unlike in [^T|] we do not consider the limit 
t — > 00. We first consider Coulomb gauge and later the general flow gauge (|3.5|). 

At lowest non-trivial order the contributions to W are due to propagators connecting 
the loop with itself. As in we ignore contributions connecting one edge with itself, 



these depend either only on t or only on R. In Coulomb gauge there are no diagrams 
which connect a space-like edge with a time-like one. Thus at second order in the fields 
there are two contributions, 

W^ 2 \t, R) = w£L(t, R) + WgL(t, R), (7.2) 
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withQ 

Wgl e (t,R) = g 2 C R [ dh [ t dt 2 (A°(t 1 ,0)A (t 2 ,m)), (7.3) 

Jo Jo 



W s ( p 2 L(t,i?) = g 2 C R [ R d Zl [ R dz 2 {A z (0, Zl z)A z (t,z 2 z)), (7.4) 

Jo Jo 

where Cr is the quadratic Casimir operator of the representation R. 
7.1 V\/ 2 ) in the leading log Langevin equation 



Calculating W(i, R) from in the theory (|1.6| ) is straightforward. The propagators are 
the ones for Eq. ( |1 . 1|) in the limit that k° ~ g A T (modulo logarithms) and |k| <C C. 
The relevant expressions are listed in Appendix |B.2| . Here 7 must be understood as the 
complete NLLO color conductivity which includes the effect of the modes with k ~ C. 
Thus we have 

A° •- • A = — 2 , (7.5) 

7p 2 

A 4 . .A' = TPl^ - - - + (p->-p). (7.6) 



obtains 



First consider W t ^ e . The integral over p is saturated by momenta |p| ~ R l , and one 



,(2) u D ^ _ C ^T± 
2vr 7#" 



W, ( ,l(i,fl) = ^^. (7.7) 



To compute W S pLe we first perform the integration over Z\ and z 2 , then the p°- 
integration which gives 

W,gL(f,J2) = -%^ f%e-^ ^ dcos^-^^-il-cosipRcose)}. (7.8) 
This integral can be simplified if we choose 

- > R 2 . (7.9) 

7 

Then the dominant contribution is from the integration region p 2 ~ 7/t < i?~ 2 . For 

Eq. (|7g ) to be reliable we need p > # 2 T, i.e., 2/7 < (g 2 T)~ 2 . Since pi? < 1 the 

11 In Ref. there is a sign error in the analogue of Eq. ( |7.3| ). Exponentiating Eq. ( |7.3| ) (with the 
correct sign) one would find that W grows exponentially for i — > oo. This error does, however, not 
affect the result for the color conductivity because the wrong sign it is used consistently. 
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Figure 2: The full A propagator in Coulomb gauge. 



square bracket can be expanded around pR = 0, and to lowest non-trivial order in this 
expansion we obtain 

Due to Eq. (|7.9| ) this is small compared to the result for the temporal contribution 

Higher loops will receive contributions from momenta of order g 2 T, and are thus 
not calculable in perturbation theory. By construction, these are the same in both 
theories (|1.1| ) and ( |1.6|) and can be discarded for a matching calculation. This is 
achieved by using dimensional to cut off both ultraviolet and infrared divergences. 
With this regularization loop corrections to Eqs. (|7.7| ) vanish. Then the leading non- 
trivial contribution to is given by Eq. ( [7.71) provided that 

g 2 T < 1 < RT l < C. (7.11) 

7.2 V\/ 2 ) in the Boltzmann-Langevin equation 

Now we compute W^ 2 -* in the theory ( |1 . 1|) and match the result with Eq. ( |7.7| ) to de- 
termine the next-to-leading log color conductivity 7. We have seen in Sec. [7[l| that 
in Coulomb gauge is determined by the A°-propagator. The lowest order contri- 
bution is again given Eq. ( |7.7|) but with 7 replaced by the 'bare' value 7b- Now loop 
corrections to this result do not vanish, and we have to compute the 1-loop contribu- 
tions to the A°-propagator. 

First let us note that the AA self-energy is zero. Every interaction vertex involves 
one A field. Thus each propagator in the loop would be a propagator which mixes A 
with a physical field. Then all propagators in the loop have singularities only in the 
lower half of the complex k° plane when k is the loop momentum. Consequently, the 
fc°-integration contour can be moved to zoo, and the integral vanishes. One can also 
understand this as a result of the invariance of the effective action T, the generating 
functional of 1-particle irreducible Green functions, under the BRST-transformation 
T is invariant because the because the transformation (12.81) is linear. A term 
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I^IP 3 ^ in T would transform into ip l U^er]^ , and there is no other term which trans- 
forms into this, so it must vanish. 

Thus the full v4°-propagator only receives contributions from the AA- and A°\- 
selfenergies (see Fig. 0). For small momenta A decouples from the / > 2 components 
of A (see Eq. ( |B.15| )), and we only have to consider the selfenergies for the I — 1 or 
vector part X 1 = C|A). The diagrams for IIaa are depicted in Appendix [B]. For loop 
momenta with k° ~ C 3 /mp, k ~ C they contribute 

TC TC 
n AA ~ g'TC- 1 x -±1 ~ log(l^)- 1 x (7.12) 

This has to be compared with the coefficient of the AA-term in the action (|4.14j) which 
is of order TC/m^. Therefore IIaa yields a contribution to the A°-propagator which 
is of next-to-leading order in log(l/gr) - . In addition to the diagrams in Appendix |D] 
there are diagrams involving the Av4°iy vertex. Using the estimates QB.5Q , (|B.12|) one 
easily sees that these are suppressed compared to Eq. (|7.12|) by powers of g. 

The contributions due to loop momenta k° ~ k ~ C are easily estimated since then 



all propagators are of the same order of magnitude (see Eq. ( B.5 )). One finds 
IIaa 



g 2 TC TC 2l . . . TC . . 

2 x — ~ 9 2 log(l/g) x — , (7.13) 



loop frequency~C 



which is suppressed compared to (|7.12| ) by two powers of g. Thus, unlike for the AA- 



selfenergy to be discussed in Sec. ||, the loop frequencies of order C can be neglected 
here. 

Now consider the AA°-selfenergy. For the diagrams listed in Appendix [D| we estimate 
n^AA ~ 9 2 T for loop momenta with k° ~ C 3 /mp. For loop momenta with k° ~ C 
we find ITa^o ~ g A T modulo log(l/g). Again the contributions which contain the 
A^PU-vertex are suppressed. 



Thus we can use the simplified form of the propagators from Sec. |B.2j to compute 
the loop integrals. For the external lines we can use the approximations in Sec. |B.2.2| . 
In Coulomb gauge we find that the A propagator for p° <C C 3 /m^ and |p| <C C is 
given by 

• = •- • + --^r x • ..|§§§ ... + 4d— 5-^-7 x 

p 4 ^IP' mf) p 4 



2T 1 AdciT p l d!i\ AQ 



+ — IIaa(O) 



7bP 2 P 2 m r> P 4 dp 



pi (7.14) 

p=0 
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o o 

Figure 3: Examples for additional contributions to the U\\ and n Aj4 o self-energies in 
the flow gauge ( |3.3|) with k / 0. 

2T +2d^|^ / ilp t ' m (k)(°|{)^G(0,k){; fe {; m | ) 



7bP 2 m^p 2 y k I k 2 



2 G t (0,k) , A d-l Cl 



Comparing this with the corresponding expression for the NLLO Langevin equation 
O we find 

- = 1 + ^1 / jlp t ^(k)(°|^G'(0,k)^ m | o ) 
7 7b mf, J k [k 2 

+ (d - l)c 2 ^^ + [1 - d Cl G t (0, k)] J. (7.15) 

for the NLLO color conductivity 7 in agreement with the result of Arnold and Yaffe 
PHI We have not checked the numerical evaluation of 7 performed in . 



7.3 Gauge fixing dependence 

In this section we address the question whether the above result for the Wilson loop 
and thus for the color conductivity depends on our gauge choice for the theory ( |1 . 1| ) . 
Consider the gauge (|3.3|) for k 7^ 0. One could leave the parametric size of k arbitrary, 
but for the purpose of power counting it is a lot more convenient to adopt a natural 
choice in terms of some scale which is already present in Eq. (|1 . 1| ) . 

In the leading log theory ( |1.6|) there would be only one natural choice. The charac- 
teristic frequency of the transverse gauge fields is k° ~ k 2 /7, which can be seen from 



Eq. ([7.6|). In Coulomb gauge the A-propagator, Eq. (|7.5|), is frequency- independent. 



12 Arnold and Yaffe ]2lJ write their results in terms of Ck = (k), which is related to our notation 
by (7k = CTrj? (k) = rn 2 3 Gt(0, k). Furthermore, their 71 is the same as our c\. 
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By choosing k ~ 7 one would make the characteristic frequencies of all gauge field 
components equal at tree level. 

For Eq. ( |1 . 1| ) there is an additional option. Consider the A propagator in Eq. (|4.31|) . 
Ge(k) has discontinuities for k° of order of the collision term C, and so does the A°- 
propagator. When k ~ C, the characteristic frequency of the transverse gauge fields 
is of order C as well. Thus one would naturally use k ~ C. This case turns out to be 
very easy to analyze. As long as k° ~ g T (modulo logs) all propagators involving A 
are approximately the same as in Coulomb gauge, and we do not have to recalculate 
any of the diagrams in Appendix [B]. Furthermore, even for k° ~ C, A is of the same 
order of magnitude as in Coulomb gauge. Therefore the suppression of the integration 
region k° ~ C found in Sec. [7.2| persists. There are additional diagrams contributing 
to the I1a a and IIaao self-energies due to the mixing of A and A, see Fig. |[ Using the 
estimates of Appendix [B] one easily finds that all these contributions are suppressed by 
powers of g. There are also contributions to the 1-loop /^-propagator of the form 

(A°A°) 1 _ loop ~ (A°A £ )U AX (XA°). (7.16) 

The longitudinal gauge field A^ is of the same order of magnitude as ^coui w hen k ~ C, 
and of order «/|k| ~ C/(g 2 T) times A° Coul when k ~ g 2 T. Therefore (A°A e ) <(A°A°). 
The self-energy IIaa is estimated in Sec. [8] as (g 2 T) 3 /m 2 :) modulo logs. Finally, from 
Eq. ( grrf ) (\A°) ~ {g 2 T)~\ so that 

(A°A)U AX (\A°) < i^l(A°A°) (7.17) 

which is strongly suppressed compared to the tree level result for (A A ). Thus we 
conclude that the result for the Wilson loop in the theory ( |1.1|) is not changed if one 
uses the flow gauge with k ~ C instead of Coulomb gauge 0. 

8 Transverse gauge field propagator at one loop 

We have seen that the NLLO color conductivity can be obtained by considering only 
the A°-propagator which gives the dominant contribution to the Wilson loop of Sec. 0. 
By simple power counting we have found that the contribution from loop momenta k 



13 Things become more complicated for k ~ 7. In this case Eq. (7.16) would contribute to the NLLO 
result for the Wilson loop. Furthermore, one would have to take into account not only self-energy 
corrections to the A°-propagator but also vertex corrections. 
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k 



Figure 4: One loop contributions to the A- A self-energy which receive equal contribu- 
tions from loop momenta with k° of order g 2 T and of order g A T. 

with [j k° ~ g 2 T is suppressed by powers of g. Thus it appears that the dynamics 
associated with this frequency scale (see Ref. fi~5fl ) is irrelevant to the soft gauge field 
dynamics. From the discussion in Sec. ^ it should be clear that one could have chosen 
another quantity for the matching calculation which receives contributions from spatial 
gauge fields as well. 

In this section we will try to see whether the spatial gauge field propagator is affected 
by loop momenta with k° ~ g 2 T, when external momenta are of order p ~ g 2 T, 
p° ~ g 4 T. It cannot be calculated perturbatively, but one can calculate the contribution 
from loop momenta large compared to g 2 T, which is the philosophy of Ref. [|TS| . 

As we discussed in Sec. [7.2| the AA self-energy vanishes. Thus we only obtain con- 
tributions which contain the AA and AA self-energies. The latter has already been 
considered in Sec. []. Consider now the contribution of the AA self-energy depicted in 
Fig. [T4] and use the estimates in Eq. ( |B.12| ). The A A propagator is of order C~ l and 
m^/C 3 for k° ~ C and k° ~ C 3 /mp, respectively. Therefore for both these frequency 
scales the product k°(AX) is of the same order of magnitude k°(AX) ~ 1. The other 
propagator in the loop is of the same order for both frequencies, (WW) ~ T/(mpC). 
The integration over 3- momenta is of order C 3 , and there is a factor g 2 from the ver- 
tices. Thus we obtain Haa ~ g 2 TC 2 /mf^. This would give a contribution to the AA 
propagator of order (cf. Eq. ( |B.18| )) 

(AA)i_ loop (p) ~ (AA>n AA (AA> ~ 4fF< AA > ( 8 ' X ) 

g 2 T 

for p ~ g 2 T, p° ~ C (g 2 T) 2 / m 2 ^ . This would mean that the 1-loop contribution is 
larger than the tree level by one power of log(l/g). However, we would expect that 

14 In this section we ignore log(l/g) in the order of magnitude estimates, so that, e.g., the collision 
term C is estimated as C ~ g 2 T. 



22 



IIaa(p) vanishes for p — * 0, and that one can expand 

n AA (p)~ P 2 n^A + ---- 

Since we consider loop momenta k with k ~ C, IT is of order 11/ C 2 . This would give 
the estimate 

(AA)i-b,(p) ~ ^-(AA). (8.2) 



In other words, the diagram in Fig. [L4] would contribute at next-to-leading logarithmic 
order. This is puzzling since the next-to-leading order color conductivity 7 which 
was calculated from the A propagator alone does not receive such contributions. It 
was argued in Sec. ||, based on the fact that Eq. (|1.1|) generates the same gauge field 
thermodynamics as Eq. ( |1.6|) , that it is sufficient to perform a matching calculation of 
one single object to obtain 7. Now it appears that it depends on the choice of this 
object whether one obtains contribution from the frequency scale k° ~ C or not. 

It is of course possible that there is a cancellation between various diagrams con- 
tributing to the AA self-energy so that the contribution discussed above come with a 
vanishing coefficient, but I do not see a reason for such a cancellation to occur. 

9 Summary and discussion 

We have found that the Boltzmann-Langevin equation (|1 . 1| ) is not renormalizable by 
power counting. There are divergent 1-loop diagrams for which the action ( |4.2j ) con- 
tains no counterterms. Introducing the necessary counterterms would invalidate the 
equivalence of the path integral (|4.1| ) and the Boltzmann-Langevin equation ( |1 . 1|) . The 
divergences appear in diagrams with A and W propagators, which for large momenta 
k fall off only like 1/k. One should emphasize that these divergences are not related to 
the ones in classical Yang-Mills or in the hard thermal loop effective theory discussed 
previously |14 |. Furthermore, they do not affect the next-to-leading log order (NLLO) 



calculations. 

We have confirmed Arnold and Yaffe's result for the NLLO color conductivity 7 
by matching the results for Wilson loops obtained from Eq. ( |1 . 1|) and Eq. ( |1.6| ) using 
Coulomb gauge. We have checked the gauge fixing independence in a certain class of 
flow gauges. The dominant contribution is due to the 1-loop A°-propagator, which 
turned out to be unaffected by loop momenta with k° 3> g 4 T[log(l/g)] 3 . 
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The gauge field propagator in the theory ( |1 . 1| ) has a non-trivial analytic structure for 



momenta of order gT, which lie outside the range of validity of (O). The calculation of 
the NLLO color conductivity is performed with an ultraviolet cutoff /x, g 2 T <C fi <C gT. 
The proof that Eq. ( |1.1| ) correctly reproduces the thermodynamics of k ~ g 2 T Yang- 
Mills fields [|15| does, however, not assume the presence of such a cutoff which could 
be a reason for being concerned. 

Finally, we have found that, somewhat surprisingly, the spatial gauge field propagator 
for momenta of order g 2 T, unlike the v4°-propagator, does receive NLLO corrections 
from loop momenta with k° ~ g 2 T\og(l/g). 

Acknowledgments. I would like to thank Poul H. Damgaard, Edwin Laermann, 
Kari Rummukainen and Daniel Zwanziger for useful discussions. This work was sup- 
ported by the DFG, grant FOR 339/2-1. 

A Properties of G and related functions 

G, defined in Eq. ( |4.15| ), is symmetric since it is the inverse of a symmetric operator, 

(v 1 |G(A;)|v 2 ) = (v 2 |G(A;)|v 1 ). (A.1) 

Furthermore, 

G{-k)=G*{k). (A.2) 

Now multiply 

[C-iv- k]G{k) = 1 (A.3) 
by (°| from the left. The collision term drops out since its Z = eigenvalue vanishes, 

-i(°|u-JfeG(ife) = ( |. (A.4) 

Now we use 

(<V = / ( V |^ = / V »(v\ 



which gives the identity 



V J V 



(1 (A.5) 



k l ( l \G{k) = -i(°\ + k°(°\G(k). (A.6) 
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From this we obtain the following relations 



and 



G £ 
k l G° l 

G 00 
k { G 0i 



A;°G 00 - 1 



1 

k 2 



;G e + i 



When k° <C |k| we can approximate G 0l (k) by 
The following matrix elements are needed to obtain the results in Appendix [D] 



f |G(0,k)|'> = p 



Cl G t (0,k)-- 
F^^'|G(0,k)|°) = Cl - ^k 2 G 00 (0,k). 



VPg^k) + k j Pl l (k) 



(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.11) 

(A.12) 
(A.13) 



B Approximated expressions and power counting 
for the propagators 

As discussed at the end of Sec. ^ we will not use the propagators of Sec. [| to compute 
the NLLO color conductivity, but instead the approximated expressions for momenta 
small compared to mr>, which are much simpler. We have to consider two cases. The 
first is k°, k ~ C. These are the characteristic momenta for the W fields which perform 
damped oscillations in such a way that the color current vanishes (see [I5|). The second 
case is k° -C C for the characteristic frequencies of the gauge fields k° ~ g A T (modulo 
logarithms of g) . There we have to distinguish k ~ C and k< C. 

B.l fc°~k~C 

When ko is of order C we have koG t ~ 1. Then one can neglect k 2 in A(k), Eq. (|4.23|) , 



mf 



i 1 
ko G t 



d 



x 



(l + OQ? 2 )) 



(B.l) 
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so that 



T 



m 



2 r t 



1 1 



+ (k -> —A;) 



(l + 0(^ 2 )). (B.2) 



Furthermore, 



and 



|5 i > = J?^G|0x(l+O(^)) 



5 ~ G - i—Po - GfG^tfti { j \G - G^)Pi 3 { j \G. 



(B.3) 



(B.4) 



To estimate loop diagrams we need order of magnitude estimates for the various prop- 
agators. From the above expressions one finds 

(A»\), (WX) ~ C- 1 
(A»A U ),(A»W),(WW) ~ Tm£C~ x ~ L^g^T- 2 {k , |k| ~ C). (B.5) 

B.2 

When the frequency k° is small compared to C the gauge field propagators are the 
same as in Ref. |2lfl . We can approximate 

1 



A ~ A 



k 2 - im 2 D k G t {0,k) 



(k ~ g*T) 



(B.6) 



\S*)=m 2 D P; j A G(0,k)\ J ) {k ~g*T) 



(B.7) 



(S°(k)\ = (°\G (k)(l-P ) 



(B.8) 



T 



2— G 00 (0,k)x (l + 0(^ 2 )) (* b ~0 4 T) 



m 



(B.9) 



(W(k,v)A (k')) 



(2ir) d+1 5{k + k')Ar(v\ (1-Po-A) G(0, k) + G(0, -k) | u ) (B.10) 



S-(fc) = G (k) - G P + PoGo + G°°P + im* VW^olW < J |Go (B.ll) 
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B.2.1 k ~ C 3 /ml and k~ C 

For these momenta the propagators are of order 

(AA) ~ mlC- 3 

(A°X),(WX) ~ C- 1 

(AA) ~ Tm 2 D C- 5 ~ L- 5 g- 8 T~ 2 

{AW) ~ TC" 3 ~ L- 3 g~ 6 T~ 2 (B.12) 

(A A°),(A°^),(^) ~ TrrtfC- 1 ~ L^g-^T' 2 . 

B.2.2 k~# 2 T 

When k is small compared to C, the propagators for the gauge fields become the 
same as in the leading log effective theory ( |1.6| ). We need the expansion of G(0, k) for 
|k| <C 

G(0,k) = ^(c l -iv-k)P (c 1 -iv-k) + (l-P )C- 1 (l-P ) 



d 



c 2 k 2 



'2 



P (v • k) 2 + (v • k) 2 P + 4l^(^ • k ) 4i? o + 0(k/C 2 ). (B.13) 

c 2 k 4 



In particular, 



and 



G(0, k)|°) = A( Cl - zv ■ k)|°) + OiC' 1 ), (B.14) 



(^WH-^CI+O^- 1 ), (B.15) 



G 00 (0,k) = ^ + O(^ 1 ), (B.16) 
k z 



G t (0) = (B.17) 



The propagators are of order 

(AA) ~ mW{g 2 T)- 2 
(A°X) ~ (g 2 T)~ l 
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1 ' C{g 2 Tf 
(A A ) ~ Lg- 4 T- 2 - TCl 



m 2 D (9 2 T) 2 

(A°w) ~ iV 4 r 2 -^- 

(HW) ~ — ^- ( Jfe ~ ^ 2 T) , k ~ p 2 T ) . (B.18) 



C Frequency integrals 



For the frequency integration it is convenient to write Ao, which is defined in Eq. QB.6Q , 

as 



where 

k 2 



(C.2) 



m 2 G t (0,k)' 

The frequency integrals for the diagrams in Appendix [D| are 

J ko y [A (k)-A (-k)] = ^, (C.3) 

l^M-k) = ^^y, (C4) 

hjslMQ - A (-k)] 2 = -ml^V. (C.5) 

We also need frequency integrals with external momentum p M = (0, p), 

/ A (k)A (p-k) = \ — — — — -i (C.6) 

A" m o k 2 G t (0, p - k) + (p - k) 2 G t (0, k) 

J o A (k)^[A (k-p)-A (p-k)] 



zG t (0,p-k) 



-, (C.7) 



A- 11 



(p - k)^ [k 2 G t (0, p - k) + (p - k) 2 G t (0, k) 
l[Ao(£-p)-A (p-*)] = ^-^. (C.f 
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D One loop results 

Here we list the Coulomb gauge results for the diagrams for the A A and A A- 
selfenergies in the limit of small external momenta p. 

One needs the expansion of Pl 3 (k — p), where k is the loop momentum, for small p, 

Vpf (k) + tiP? (k)l + 0(p 2 /k 2 ). (D.l) 



p 



P^(k- P )=Pf(k) + ^ L 

For p — ■> the A A-selfenergy becomes U l £ x (p) — > 5^U\x(0). The diagrams con- 
tributing to — IIaa(O) are 

,2 




ni 




D 

Ng 2 T 2 



mt 



2d 
k 



k 4 G t (0,k) 



2 P t im ( M |G(0,k)| fcm ) 



(D.2) 



n 2 



d(<1 11 k*G,(0,k) 




2(rf-l)^ 2 T 2 /" G 00 (0,k) 




Ng 2 T 2 



VTLt 



k 2 



D 



_L( Cl _l k2G oo (0 , k ) 




0. 



(D.3) 

(D.4) 
(D.5) 
(D.6) 
(D.7) 



Note that the diagram ( p.6| ) enters IT AA with an additional factor 2 because the reflected 
diagram gives the same contribution. 



Arnold and Yaffe [211 write their results in terms of a different matrix element than 



the one Eq. ( p.3|) . Ours is related to theirs by 



Pl m { kl \G{0,k)\ km ) = P! m (°\v k v l G{0,k)v k v m \°) + , 

ak 2 



2ci d + 1 
d 2 



G oo (0,k). (D.8) 



With this relation the sum of Eqs. ( |U.2| ) ( [DTB| ) gives 
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n AA (0) = 2d^Xf^ I s - 
ml .L k 2 



P^(k)(O\vWG(0,k)v k v m \°) + (d- l)c 2 - Gt(0 ' k) 



k 1 



(D.9) 



The Av4°-selfenergy Tl\ A0 (p) vanishes for p = 0. Expanding around p = the lowest 
order contribution in p is given by the diagrams 

../~V = NfT^f' (D.10) 



k 




+ 



a- 




-Ng 2 T(d - iy J i ( Cl G t (0, k) - ~) (D.ll) 
-Ng 2 T Cl (d - iy jf ^G t (0, k) (D.12) 



the sum of which gives 



n^o = ■ • • ^jjk - = 2Ng 2 T(d - iy J 1 Q - Cl G t (0, k)^ . (D.13) 
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